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Abstract 


A nonlinear circuit model of a polymer electrolyte membrane (PEM) fuel cell stack is presented. The model allows the simulation of both 
steady-state and dynamic behaviour of the stack on condition that the values of some of its parameters are changed in the two operating conditions. 
The circuit parameters can be obtained by means of simple experimental tests and calculations. A commercial PEM fuel cell stack is modelled as 
seen from the power conditioning system side, without requiring parameters necessary for complex mathematical models and not easily obtainable 
by the majority of users. A procedure of parameter determination is developed and a comparison between the simulated and experimental results 


for both steady-state and dynamic behaviour of the PEM stack is shown. 


© 2006 Elsevier B.V. All rights reserved. 
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1. Introduction 


Rising crude oil prices (with highs over US$ 70 a barrel dur- 
ing the last year) pose the exploitation of alternative energy 
sources as a serious challenge for the next future. Fuel cells 
appear one of the most appealing renewable energy technolo- 
gies for their low environmental impact and high conversion 
efficiency. Their application ranges from stationary to portable 
power generation, including transportation. Modelling fuel cells 
is then necessary to simulate the behaviour of more complex 
systems (e.g., electric vehicles, or electric low-power plants 
or cogeneration systems), in which fuel cells are integrated as 
source of energy. Lots of papers present mathematical models 
for PEM fuel cells, in which their typical application require- 
ments of high specific power, rapid start-up, low-temperature 
operation and ease of construction are met. Different load con- 
ditions, temperature and pressure of gases, as well as spatial 
dimensions of the cell [1], can thus be taken into account at 
the design level and simulation results can be of help in setting 
up operational strategies. Most of these models are however 


* Corresponding author. Tel.: +39 051 2093484; fax: +39 051 2093588. 
E-mail addresses: ugo.reggiani ? mail.ing.unibo.it 
(U. Reggiani), leonardo.sandrolini @ mail.ing.unibo.it (L. Sandrolini), 
gianlorenzo.giuliattini € mail.ing.unibo.it (G.L. Giuliattini Burbui). 


0378-7753/$ — see front matter © 2006 Elsevier B.V. All rights reserved. 
doi:10.1016/j.jpowsour.2006.11.062 


extremely complex, involving many partial differential equa- 
tions and their boundary conditions, and need a lot of expertise 
in identifying and estimating their large number of parameters 
as detailed information on the cell is essential [2-4]. The variety 
of material properties to know, such as porosity, permeability, 
effective diffusion and charge transfer coefficients, makes these 
models cumbersome, and not often easily exploitable. Simplified 
models have also been proposed, in which the reversible volt- 
age and the voltage drops are summarized in simpler equations 
[5-7]. On the other hand, complex mathematical models can 
be simplified with the introduction of some empirical equations 
instead of partial differential equations [8,9]. These equations 
use fitting coefficients obtained from experimental data and are 
therefore related to a particular operating condition. This means 
that the equations may fail to predict experimental data in dif- 
ferent conditions. The steady-state performance of a fuel cell 
stack can be represented by the majority of these models, result- 
ing in the so called polarisation curve, that is a plot of voltage 
versus current density for a given set of operating conditions. 
Besides, the dynamic behaviour of a fuel cell stack cannot be 
disregarded in all those applications where mechanical, thermal 
or electrical quantities have fast variations [9,10]. Commercial 
fuel cell systems are power modules that include also auxil- 
iary circuitry to control subsystems for fuel and air supply, and 
water disposal. This paper will be focused to model a commer- 
cial fuel cell stack as seen from the power conditioning system 
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Fig. 1. Fuel cell stack used in the experiments. 


side, i.e., electronic converters which interface such source to 
the electrical grid and/or the load, without requiring param- 
eters either not easily obtainable by the majority of users or 
confidential to the manufacturers only. The model represents 
the overall system, i.e., fuel cell stack and auxiliary circuitry, 
as shown in Fig. 1. The performance of the stack is suitably 
described with a simple circuit model that can be represented in 
terms of an equivalent circuit. The model allows the simulation 
of both steady-state and dynamic behaviour of the stack. With 
reference to other circuit models of PEM fuel cells that have 
been recently presented in the literature [11—15], the parame- 
ters of the proposed model can be determined through simple 
experiments and calculations which still offer insight into the 
nonlinear phenomena occurring in a fuel cell. Moreover, the 
model is applicable to PEM fuel cells that lack detailed informa- 
tion. A comparison between experimental and numerical results 
obtained with the model presented in this paper is made and 
discussed. 


2. Basic operating principles of a PEM fuel cell 


A PEM fuel cell consists of two electrodes separated by a 
thin layer of catalyst in contact with a plastic permeable mem- 
brane that allows protons to pass through but prohibits the 
passage of reactant gases. The electrodes and the membrane are 
pressed between two conductive plates containing some chan- 
nels through which the gases (hydrogen and oxygen) are fed 


from gas supply chambers. The electrodes and the membrane, 
which acts as an electrolyte, form a structure with an overall 
thickness of less than a millimeter. The electrodes have a porous 
structure which allows the reactants to flow. The chemical reac- 
tions that occur at the catalyst layer yield hydrogen ions and 
electrons at the anode 


Hz 2H*-42e', (1) 
and water (and heat energy as a by-product) at the cathode 
10» + 2H* + 2e7 > H20. (2) 


The hydrogen ions cross the membrane, while the electrons 
flow through the external circuit connected to the electrodes gen- 
erating useful current. To achieve the power and voltage required 
by the system connected to the fuel cell, several single cells can 
be connected in series into a fuel cell stack. 


3. Circuit models 


Irreversible losses that occur under operating conditions 
make the voltage of a fuel cell stack less than its reversible value, 
Viey. This is the theoretical value of the open-circuit voltage, and 
for a hydrogen fuel cell stack it is given by the formula [5] 

Agr 

Vrey = —Neell OF’ (3) 
where gf is the molar Gibbs free energy, F the Faraday con- 
stant and Nee is the stack fuel cell number. A simple analytical 
model can be based on the main voltage drops only, i.e., the 
activation voltage drop Vact, the ohmic voltage drop Vohm, and 
the concentration voltage drop Veonc. Voltage drops in a fuel 
cell can also occur either for cathode flooding, when the water 
produced by the chemical reaction is not effectively removed 
[10], or for membrane dehydration. In many fuel cell models 
the membrane is considered as being completely saturated of 
water [9] and therefore hydration issues of the membrane are 
neglected. Besides, when all the variables that govern the rate 
of water formation (i.e., the current) or removal (i.e., the tem- 
perature and pressure) are kept in the range prescribed by the 
cell manufacturer, also the voltage drops due to cathode flooding 
may be ignored. Other causes of voltage drop are related to the 
diffusion of some hydrogen molecules and of some electrons to 
the cathode through the membrane. The former phenomenon, 
known as fuel crossover, is more important, and represents a 
waste of fuel. Both phenomena are basically equivalent, as they 
result in two electrons per hydrogen molecule that go from the 
anode to cathode internally rather than through the external cir- 
cuit, and can be referred to as an internal current. Although the 
amount of electrons is small, it does cause a sensible voltage drop 
at open-circuit, especially with low-temperature cells, such as 
PEM ones [5]. As this voltage drop is particularly significant 
for small currents, it can be comprised in the activation voltage 
drop. The total stack voltage, Vr, is then obtained excluding all 
the voltage drops from the reversible voltage 


Vr = Vrev m Vact = Vohm aud Vconc- (4) 
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Ract 


Cg 
Fig. 2. Equivalent circuit of a PEM fuel cell stack. 


The dynamic behaviour of a PEM fuel cell stack can be rep- 
resented by means of a simple first-order equivalent circuit such 
as the one depicted in Fig. 2 [5,7,16]. The ohmic voltage drop 
is represented through the resistor Rohm, which expresses the 
internal resistance of the cell stack, i.e., the resistance due to 
nonideal electrodes and conductive plates and to proton transfer 
through the membrane [5,9]. The activation voltage drop is rep- 
resented through the parallel connection of a resistor, Ract, with a 
capacitor, Ca, that models the double layer of charge at the inter- 
faces between the membrane and the electrodes. The dc voltage 
source Vo is the effective open-circuit voltage of the fuel cell 
stack. The model does not provide for the concentration volt- 
age drop, related to changes in the concentration of oxygen and 
hydrogen at the electrodes due to the continuous chemical reac- 
tion. Although in a well designed system with good oxygen and 
hydrogen supplies this drop should be very small at the rated cur- 
rent [5], it may be taken into account by an impedance in series 
with Ract. As the elements of this simple circuit are assumed lin- 
ear, the simulation of complex nonlinear phenomena such as the 
ones occurring in a fuel cell stack is not always appropriate. In 
particular, neither the internal current effect, nor currents greater 
than the rated one can be treated with this model. Furthermore, at 
steady-state the equivalent circuit of Fig. 2 is reduced to a resis- 
tive circuit and, consequently, the voltage drops at steady-state 
are simply proportional to the current supplied. The voltage ver- 
sus current density curve of the fuel cell stack would therefore 
result in a linear steady-state performance curve, that approxi- 
mates the experimental performance curve just for intermediate 
current density values, but not for low or high current density 
values. The aim of this paper is to present an equivalent circuit 
able to approximate both dynamic and steady-state behaviour of 
the stack. 


Fig. 3. Proposed equivalent circuit. 


4. Proposed model 


In order to simulate the cell stack behaviour effectively, in 
particular in steady-state conditions, a more complex circuit 
model is developed, with nonlinear active and passive compo- 
nents. With reference to Fig. 3, the circuit contains two nonlinear 
current controlled voltage sources Vcsi(/) and Vces2(/), which 
can be related to the activation and concentration voltage drops, 
respectively. Asitis known, these voltage drops can be expressed 
as 


I+ ihn 
Vact = fi (T) = Astack In ( n ) , (5) 
0 
and 
I 
Vconc = f20) = —B In ( ni ) . (6) 
I 


In (5) Astack is Tafel ’s stack coefficient; Z, In and Jp are the output, 
internal and exchange currents of the stack, respectively. In (6) 
B is an empirical coefficient and A is the limiting current. It can 
be assumed that Vcs2 2 f2(I). 

With respect to the model proposed in [14], the developed 
model considers the effect of the internal current on the acti- 
vation drop and takes the influence of the temperature on the 
reversible voltage into account [2,6], being 


Viey = Neel E — 0.9 x 1073(T — 298) 


1/2 
RT 2 

i| ee (7) 
2F PHO 


where R is the universal gas constant, T the operating absolute 
temperature and p is the partial pressure of the different species. 
Constant pressures of gases are assumed to simplify the model. 

In regard to the passive elements of the equivalent circuit 
shown in Fig. 3, the capacitor C affects the transient behaviour 
whereas the meaning of the resistors R; and R? is highlighted in 
the next section. 


5. Experimental 


The analysis of the steady-state and dynamic behaviour of the 
model is made by comparing simulated and measured results. 
The experimental tests were carried out on a commercial PEM 
fuel cell system composed of a cell stack and auxiliary circuitry 
for internal regulation. The stack was a series of 47 membrane 
electrode assemblies (MEAs) comprising a Nafion 112 mem- 
brane with a thickness of 50 jum. The active area of the electrodes 
was 115.8 cm?. The system had a rated output power of 1.2 kW 
at a voltage of 26 V and the open-circuit voltage was 43 V. The 
auxiliary circuitry shuts down the fuel cell stack when the stack 
current reaches the rated current value. 

The electrochemical impedance spectroscopy (EIS) is a 
widely accepted method for parameter determination of PEM 
fuel cell stack equivalent circuits. In this tests small signals 
of variable frequency alternating current are superimposed to 
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a steady-state operating point of the cell stack. The fuel cell 
stack impedance spectrum is obtained from the measured volt- 
age across the stack, and the parameter values of the adopted 
equivalent circuit can be determined using a fitting procedure. 
Nevertheless, this method cannot be applied to equivalent cir- 
cuits with active elements, such as controlled sources. In this 
case, the current interrupt method (CI) is well suitable [5,16—19]. 
When the current of the cell stack is interrupted, the stack volt- 
age has an instantaneous rise, A Vonm, equal to the ohmic voltage 
drop, and then it moves to the open-circuit voltage value, Vo, with 
an increase, A Vac, that equals the change of the activation volt- 
age drop, provided that the operating current is low enough so 
that the concentration voltage drop is negligible. However, it may 
often result very difficult to accurately carry out the graphical 
estimates of the voltage rise, which make the method so simple in 
principle. In fact, the point where the vertical transition ends may 
not be exactly discriminated, yielding a possible overestimate 
of the ohmic voltage rise. Another difficulty in interpreting the 
results is introduced by real oscilloscopes, which do not show the 
vertical voltage rise ideally expected. However, these limitations 
can be overcome by extrapolating the transient voltage curve at 
the moment of the current interruption with an appropriate time 
scale [16,20—22]. In this way, the ohmic contribution, AVonm, 
to the voltage rise of the stack can be separated from that related 
to the activation phenomena, AV. Several CI tests were per- 
formed at different operating conditions, which can be varied by 
connecting the fuel cell to a home-built electronic load, depicted 
in Fig. 4. Eight modules, each composed of power resistors and 
a N-type MOSFET transistor, dissipating a maximum power of 
150 W, are connected in parallel. A time range of 100 ws was 
chosen to separate the ohmic voltage rise contribution from the 
activation one (see Fig. 5) through the extrapolation procedure. 
Although voltage spikes preceding and following the interrup- 
tion, due to the feedback control of the electronic load, can be 
noticed, they do not affect the estimate of the ohmic voltage rise. 

Voltage transients following the instantaneous ohmic voltage 
rise during CI tests were measured. The whole test set-up, com- 
prising the PEM fuel cell stack and the electronic load, is shown 
in Fig. 6. 


Fig. 4. Home-built electronic load used in the CI tests. 


36 E T 1 1 


stack voltage [V] 


time [us] 


Fig. 5. Voltage transient (solid line) and fitted voltage (dashed line) following a 
10.04 A current interruption. 


Fig. 6. Test set-up used in the CI tests. 


The experimental tests were completed with the measurement 
of the steady-state stack voltage versus current characteristic. 
The stack temperature varied between 25 and 65 °C. 


6. Results and discussion 


The results of the CI tests are summarized in Table 1. For each 
value of the operating current / the ohmic resistance is evaluated 


Table 1 

PEMFC CI test results 

Test I(A) A Vohm (V) Rohm (Q) 
1 4.4 0.14 0.031 

2 6.2 0.25 0.040 

3 8.2 0.32 0.040 

4 10.0 0.43 0.043 

5 12.2 0.49 0.040 

6 14.3 0.58 0.041 

7 16.2 0.66 0.041 
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from the determined ohmic voltage rise, AVohm, as 


A V, 
Rohm = a (8) 


The value of the ohmic resistance of the PEM fuel cell sys- 
tem varied from 31 to 43 mQ in the tests. An average value 
of 39 mQ (standard deviation o=4mQ) was then chosen for 
the resistor Rohm. This value corresponds to an areal resis- 
tance of 0.096 Q cm? and a resistivity in the normal direction of 
19.22 Q cm which are consistent with published data [23—27]. 
However, it is difficult to make comparisons as the measure- 
ment techniques are different and the tests are carried out on 
a single cell or an isolated membrane in some cases and on a 
whole stack in others. With respect to [28], the lower resistance 
value found in this paper can be related to the shorter overshoot 
produced by the improved load. The absence of inductive oscil- 
lations allowed the ohmic voltage rise to be determined with 
higher accuracy. 

Recalling (5), the activation voltage rise AVac that occurs 
during the CI transient can be written as 


AVaa = fi) — fi (0). (9) 


It was chosen not to estimate A Vac from a transient output 
voltage as the remaining contribution to the voltage rise. In fact, 
a steady-state condition could not be necessarily reached by the 
fuel cell module even with a longer transient time range. The 
quantity AVac is evaluated through the experimental steady- 
state stack voltage versus current characteristic. In fact, as at 
low currents the activation loss is almost entirely responsible 
for the stack voltage drop [2.6], the initial trait of the experi- 
mental steady-state stack voltage versus current characteristic 
can be used to approximate the activation voltage drop fi(), 
provided that the ohmic voltage drop is taken into account. 
These fi(J) experimental values allow the activation voltage 
drop versus current curve to be built in the operating current 
range of the PEM stack. The least-squares method used to fit 
(5) to these experimental values yielded Agtack =3.0+ 0.1 V, 
I, 20.6 0.1 A, Jo 25 -E 1 mA. These values correspond to the 
Tafel’s cell coefficient Ace 2 63.8 mV and to the current den- 
sities Jn 2 5.2mA cm ? and Jo 2 0.04 mA cm ?, which are in 
the range of values reported in literature [5,29—32]. The acti- 
vation voltage drop fi(J/) is shown in Fig. 7. From Fig. 7 
A Vact — fi) — f1(0) is obtained. 

The dynamic behaviour of the proposed equivalent circuit is 
examined when disconnecting a load. With reference to Fig. 8, 
the time constant t of the R; C parallel can be evaluated through 
a fitting procedure of the voltage transient following the instan- 
taneous ohmic voltage rise during a CI test with a time range 
of 5s. At the moment of the load current interruption, the volt- 
age across the RıC parallel is equal to AVact. In fact, at the 
end of the voltage transient following the instantaneous ohmic 
voltage rise, the stack voltage rise is equal to AV,,; and in the 
proposed equivalent circuit this voltage rise is due to the capac- 
itor discharge through R4. As the voltage across a capacitor is 
continuous in time, the voltage across Rj is AVact also at the 
moment immediately preceding the load current interruption, 
when the steady-state current / flows through R;. Consequently, 
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Fig. 7. Activation voltage drop vs. current curve for the PEM fuel cell stack. 


R, can be calculated as 


A Vact = fi) - fi) 
I I ` 


Being t = R,C, the capacitance C can be obtained at the dif- 
ferent operating currents. Fig. 9 shows the capacitance values 
obtained with the voltage transient fitting procedure and the 
analytical fitting curve 


Ri) = 


(10) 


CU) =a In (=) +b, (11) 


ref 


where Ief is a unitary reference current, a=0.034 + 0.001 mF 
and b=0.040 + 0.003 mF. During the transient following the 
instantaneous ohmic voltage rise (4) reduces to 


Vr = Vrev E Vact (12) 


being Vohm =9. Moreover, as the concentration voltage drop is 
usually negligible for currents lower than the rated one and it 
is difficult to take it into account dynamically, the current con- 
trolled voltage source Vcs is neglected. In the equivalent circuit 
of Fig. 3 this transient is represented by the capacitor discharge 
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Fig. 8. Measured and simulated stack voltage transient following a 10.04 A 
current interruption. 
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Fig. 9. Capacitance values for the proposed equivalent circuit. 


through the resistor A1 and the Kirchhoff's voltage law yields to 


Vr = Vrev — Vesi — Ve, (13) 
where 
Vc = [fA — fi] e "6, (14) 


From (12)-(14) it follows that during this transient the fol- 
lowing equation holds 


Vac = Vesi + [f(D — fi(0)]e "€, (15) 


It is easy to verify that the value of the current controlled 
voltage source Vcs; to be considered in order to satisfy the 
boundary conditions of the transient is Vcs; —f1(0). In fact, the 
value of Vact at the moment immediately preceding the interrup- 
tion of the load current (time t=0) is then fı (D), and the value 
of Vact at the end of the capacitor discharge (t — oo, no-load 
steady-state behaviour) is f1(0). It can be noticed that at this 
last operating condition Vr = Vo and thus from (5) and (12) it 
follows 


fi(0) = Asa In (2) = Ves — Vo. (16) 

As Fig. 8 shows, a good agreement between the measured 
stack voltage transient and that simulated with the proposed 
circuit is obtained. For transients different than the one exam- 
ined, where the load current varies suddenly from J; to h 
steady-state values, the above considerations still hold, provided 
that currents 7 and O are replaced by 7, and Jy, respectively. 
It follows that for a given transient the resistor R depends 
on both bound state currents of the transient. The simulation 
for a stack voltage transient in which the current decreases 
from 6.10 to 2.39 A 1s plotted with the experimental result in 
Fig. 10. The results are still in good agreement also for larger 
current variations, as shown in Fig. 11, where the transient 
following the current variation between 19.99 and 12.03 A is 
depicted. 

At a steady-state operating condition the voltage across the 
capacitor (14) is equal to zero but the capacitor behaves as an 
open-circuit. These two conditions can be satisfied simultane- 
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Fig. 10. Measured and simulated stack voltage transient for a current variation 
from 6.10 to 2.39 A. 


ously only if at steady-state the value of R, equals zero. From 
the circuit of Fig. 3 we have 


Vr = Vev — Vesi — Rol, (17) 


where / denotes a steady-state current. Introducing the expres- 
sions of Vact and Vohm into (4) we can write 


Vp- Viev m A) — Ronmf. (18) 


Eqs. (18) and (17) yield R? = Rohm and the steady-state value 
Vcsi —2f1(). In Fig. 12, good agreement can be noticed between 
the measured and calculated steady-state stack voltage versus 
current characteristics. The concentration voltage drop was not 
taken into account in the calculation, as the limiting current J, 
was not known. This drop occurs at very high current densities, 
which could not be anyway reached by the fuel cell stack under 
test, as the maximum operating current value was purposely 
limited by the constructor. It can be noticed that the controlled 
source Vcs; gives the activation voltage drop Vact — fi (J), being 
I the current at a steady-state operating condition or the steady- 
state current at the end of a transient. The resistor R; has a 
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Fig. 11. Measured and simulated stack voltage transient for a current variation 
from 19.99 to 12.03 A. 
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Fig. 12. Measured stack voltage and voltage drops implemented in the model 
vs. current in steady-state conditions. 


resistance equal to zero at steady-state and is a nonlinear param- 
eter depending on the current during a transient. Consequently, 
the capacitor C connected in parallel with R; is also a nonlinear 
current-dependent parameter. The dependence of the resistance 
Rohm = R2 on the current may be neglected choosing its aver- 
age value in the considered operating current range of the PEM 
stack. 

It has to be highlighted that the aim of this paper is mainly 
to present a model for a commercial PEM cell stack supply- 
ing the electrical grid and/or a load through the interfacing with 
electronic power converters. The objective is to have at one's 
disposal a model for the whole complex power electrical sys- 
tem of which the fuel cell stack is part. As a consequence, the 
proposed model has been validated in the operating range of a 
commercial fuel cell stack whose auxiliary circuitry does not 
allow currents higher than the rated one to be obtained. How- 
ever, it can reasonably be expected that the proposed model is 
suitable to take account of the concentration voltage drops, too. 
In fact, the stack voltage versus current characteristic could be 
experimentally measured in a stack where the auxiliary circuitry 
can be deactivated. If the limiting current is known from the cell 
manufacturer, the model with the controlled source Vcs» can be 
fitted to the experimental results in order to find the empirical 
coefficient B. The comparison between measured and calculated 
voltage versus current characteristics allows the model to be 
validated. 


7. Conclusions 


A circuit model for a PEM fuel cell stack is presented. The 
model is nonlinear and can be used to simulate the steady- 
state and dynamic behaviour of the PEM stack if the values 
of some of its parameters are changed in the two operat- 
ing conditions. The circuit parameter values are determined 
through simple experimental tests and calculations. A good 
agreement between the simulated and experimental results for 
both steady-state and dynamic behaviour of the PEM stack is 
obtained. 
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